Density functional derived from Feynman 

diagrams. 

Georges Ripka 

ECT*, Villa Tambosi, 1-38050 Villazano (Trento), ITALY 
Institut de Physique Theorique, 
Centre d'Etudes de Saclay, F-91191 
Gif-sur-Yvette Cedex, FRANCE ripka@cea.fr 

December 22, 2009 

Abstract 

We construct a stationary density functional for the partition function 
from a chosen set of one (boson) line irreducible Feynman diagrams. The 
construction does not proceed by the inversion of a Legendre transform. 
It is formulated for fermions with Coulomb interactions, in which case a 
stationary functional of the particle density is contructed, as well as for 
nucleons interacting with mesons, which involves a stationary functional 
of several densities. The use of Kohn-Sham orbits is shown to be an 
unnecessary complication. 

1 Introduction. 

The advantage of constructing density functionals of finite systems, such as nu- 
clei, atoms and crystals, resides in the fact that they give rise to single particle 
states which can be calculated with a static (energy independent) potential. 
This makes computations considerably easier if one needs to go beyond the 
mean field approximation. For systems which are composed, for example, of 
electrons with a Coulomb interaction, it is natural to construct a functional of 
the particle density n (r) — (ip> (r)ip (r)) because this is the density to which 
the Coulomb field couples. It is for such systems that density functionals were 
first constructed [U [5] . In nuclei however, there is no reason to single out the par- 
ticle density n (r). For example, pions couple to the density (ip* (r) Taj^aip (r)), 
scalar and vector mesons couple to densities Up* (f) 7o'0 (f)) and (r) ip (r)) 
respectively. In this work, we show how to express the partition function of the 
finite system in terms of a stationary functional of densities which are deter- 
mined by the nature of the interactions. 

The construction of stationary functionals of the density matrix, appropriate 
to systems with two-body interactions, was developed in Ref.[3], based on an 



1 



earlier work 0] which described the theory of a static potential in the presence 
of two-particle two- hole excitations. In these works, the density functional 
are derived from a chosen set of one (interaction) line irreducible Feynman or 
Goldstone diagrams. More recently, considerable work has been devoted to 
the construction of stationary density functionals obtained from a Legendre 
transform (see [5], [5], [7j,[S], 0,110] and references therein). 

In this work, we construct a density functional from a chosen subset of 
one line irreducible Feynman diagrams which are assumed to be dominant. We 
consider both fcrmions with Coulomb interactions and nucleons interacting with 
scalar, vector and 7r mesons. The formalism is simpler in the former case because 
it involves only one single-particle density, whereas several densities are involved 
in the latter. The partition functions of the two systems are defined in Sections 
[2] and [3] The Feynman diagram rules for the two systems are summarized in 
Appendices [A] and [BJ Section |4] describes the properties of one-line irreducible 
diagrams. In section [5] the single particle densities are expressed in terms of 
one-line irreducible diagrams which are calculated in terms of self-consistent 
potentials. The potentials and single particle densities are related by a simple 
analytic expression. Section [6] describes iteration procedures to calculate the 
self-consistent potentials and the single particle densities. In Section [7| the 
partition function is expressed in terms of one-line irreducible diagrams and in 
Section [5] it is shown to be a stationary functional of either the potentials or 
of the single particle densities. In Section [TU] the theory is expressed in terms 
of so-called Kohn-Sham orbits the use of which is shown to be an unnecessary 
complication. 



2 Non-relativistic fermions with Coulomb inter- 
actions. 

Non-relativistic fermions, such as electrons with Coulomb interactions, can be 
described by a hamiltonian of the form: 

H-nN= Jd 3 r (f) (~ + U ext (r) (r) 

+ \ [ d 3 n d 3 r 2 ^ {fx) (r 2 ) 6 ^ (r 2 ) ^ (fi) (2.1) 

where if) is the fermion field, M the fermion mass, e its electric charge and 
Uext (r) an external potential which could be due, for example, to an applied 
external field, to a nucleus or to the crystal nuclei. We define: 

ho = + U ext (r) (2.2) 

The partition function is given by: 

Z = Tr e -«*-/dV) = e w = f D ^ e -r(^>) (2 .3) 



2 



where tp^ and tp are independent Grassmann variables and where I (ip^,ip\ is 
the euclidean action: 

ifi (a;) (d T + h - (i) ip (x) — - [ d i x l d i x 2 tft (x x ) ip (x x ) {x\ \K\ x 2 ) V' t (^2) ip (x 2 ) 



" 2 J ' ' ' -------- ■ - 

(2.4) 

with x^ = x^ — (r, f) and: 

-K=4j2 -(^W^)^(n-r 2 ) 4H ^ f2| (2.5) 

Note that we defined the Coulomb potential to be — K (and not K) so that 
K is a negative definite operator. The reason is to use a common notation 
for fermions with Coulomb interactions and nucleons interacting with mesons, 
described in Section [3] We can also express the partition function in terms of 
an integral over an additional boson field u)q- 

Tr e -KH-*N) = e w = J D(^,i>)D (coo) e -'(^-**0 

I (ip\ip,uj ) = J d 4 x ipl (x) (h — fi + iu>o) J d A xid A x 2 uj (xi) (xi l^" 1 ] x 2 ) uj (x 2 ) 

(2.6) 

By integrating over the field luo we recover the action (|2.4[) . By integrating over 
the nucleon fields ip* and ip we obtain the partition function in the form: 

Z = Tre-K H -W =e w = ( D (lu ) e' 1 ^ 



I (ujo) = -Tr In (d T + ho - fi + iu )—^ J d 4 x 1 d i x 2 uj (xi) (x x \K 1 \x 2 )wo(x 2 ) 

(2.7) 

In (|2.6p the trace is taken over the Hilbert space composed of 0, 1, 2, ... fermions. 
In equation (|2.7[) and in the following, Tr is a trace over the euclidean space-time 
(r, r) and over the spins of the fermions, in other words, over the arguments of 
the fermion field ip: 

Tr A = tr J d 4 x (x \A\ x) = dr J d 3 r ^ (rfa \A\ rrer) (2.8) 

and tr denotes a trace over all the discrete indices other than (r, r) : 

trA = Y,{°\A\v) (2.9) 

(T 

The chemical potential of the fermions is \i. In the zero temperature limit 
pi — > 00, the energy E of the system composed of N fermions is given by: 

E — fiN — ~t;W (2.10) 

P 

The reader who is not interested in nucleons interacting with mesons can 
proceed to Section 2] 
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3 Nucleons interacting with mesons. 

We consider a system of nucleons of mass M interacting with a scalar meson a 
of mass m a , an isoscalar vector meson lo^ of mass m u and pions 7r a of mass m w . 
The partition function can be expressed as the following path integral over the 
nucleon and meson fields: 



Tre -/3(H-„N) = £ w / D ( ivtj iv) D (a) D («) D (tt) g-^W.".*) (3.1) 

where l(N f ,N, a, tt, u) is the euclidean action: 

Ij {N\N,cj,it,u) 

= iV f ^d T + + /3M -fi + g(3a + ^j-^Taj^ {d^ix a ) + jp0s a bcr c 7ra"/n (d^b) + fl^/H^l N 

/\ 1 1 1 

d A x +-a (-d 2 + to 2 ) <t + -vr Q (-<9 2 + to 2 ) TT a + - W/4 ((-<9 2 + m 2 ) ^ + fyft,) w> 

(3.2) 

The euclidean action (|3.2[) is expressed in terms of the following euclidean 
(fi^ = <V) 4- vectors: 



c fl =x^ = (r, r) 7m = 7 M = (*A f) w m = ^ = ( w o, 0) 
a 2 = a 2 + V 2 y d 4 x = jf " dr J d 3 r 



(3.3) 



We shall work with the euclidean action obtained by integrating out the nucleon 
fields, in which case the partition function becomes: 



(3.4) 



with: 



a ■ V 



J(ct,7T,w) 



fjA 



I 



Tr In <9 T H ; 1- (3M - fi + g(3a + — — /?T a 7 5 7 M (<9 M 7r a ) + — ^/fea&cTcTT^ (<9 M 7r fc ) + g^P^^uj 



4/ 2 ' 



-C7 (-<9 2 + m 2 ) a + -tt (-<9 2 + to 2 ) vr a + -0J„ ((-d 2 + to 2 ) ^ + d^dj) . 

(3-5) 

The trace Tr is taken over the quantum numbers which define the fermion field 
N (x). For example: 



Tr ( d T + + (3M - n + g(3a J = / d A x tr ( x 



a ■ V 

d T H h /3Af - /i + g^cr 

z 

(3.6) 
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and tr is a trace over the discrete (Dirac, flavor) quantum numbers other than 
X = (r, r). 

We stress from the outset that we are not suggesting that (|3.5p is the action 
which should be used to calculate finite nuclei. We are simply introducing the 
a, lo and 7r mesons to show how to construct a density functional with different 
types of meson-nucleon interactions. 

In order to express the partition function in terms of Feynman diagrams, we 
adopt a condensed notation. We define the unperturbed fermion hamiltonian: 

h a = ^-^ + [3M (3.7) 

i 

The following is not restricted to a Dirac unperturbed hamiltonian. It could, 
for example, be replaced by its non-relativistic reduction ho = —^jj , in which 
case the couplings would have to be modified accordingly. The hamiltonian ho 
acts in the Hilbert space of one fermion. We define the following T a operators 
which act in the same space: 

r a = (rw rH,r«) = (g^g^, (3.8) 



®a(U = tWW-c7i. (3-9) 



as well as the operator: 

1 

The meson fields are denoted by: 

S a (x) = (a (x) , (d^ a (*)) , ^ (x)) (3.10) 
We will use the condensed notation: 

(rs) = rW<7 + rM W(t + rW(Va) 

= g a 0a + ^-/3r a 757 Al {d^ a ) + g^h^i^ 



S (95) = Tr a <3>a(pb) (dfiTTf,) = jj^f3E abc T c Tr a 'y IJ , (<9 M 7r fc ) (3.11) 



1 

In the expressions above, 7^ and lu^ are the euclidean 4-vectors defined in (|3.3 
Note that the coupling constants are included in T and 0. 
We further define the meson propagators K: 

1 



-d' A +ml 

K=\ 5 ab _ & } +mi I (3 12 ) 

With this condensed notation we can write: 

^SK-'S = f d A x j X -g {-d 2 + m 2 a )a + ^ir a (-d 2 + ml) n a 
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+5 w m + m " ) <W + d„du) w v | (3.13) 

and the euclidean actions (|3.2p and (|3.5[) acquire the form: 

I (A rt , AT, 5) ee ATt («9 r + h — fi+ {TS) + 5 (95)) JV + ^SK^S 

7(5) ee -Tr In (9 T + /i -M + (r5) + 5(95)) + -SJf-^ (3.14) 
The partition function is given by: 

Z = Tre^ l3{H ' tlN} ee e w = [ D (n\ N) D (5) e -NHa r +h -^+(rS)+S(es))N-iSK- 1 s 



= j D{S)e TrXa{dT+ha '^ +{rs)+s[es)) ~^ SK ' ls (3.15) 

The inclusion of separate chemical potentials for neutrons and protons would 
cause no difficulty. In the zero temperature limit (3 — > oo, the energy E of the 
system composed of A nucleons is equal to: 

E-fj,A = — -W (3.16) 

P 

4 One (boson) line irreducible diagrams. 

In this section we define the one (boson) line irreducible diagrams which are 
used to construct the density functionals. The rules for calculating Feynman 
diagrams are summarized in Appendices \X\ and [Bl Consider a general unlabeled 
connected diagram, such as: 





(4.1) 

We can recognize so-called articulation lines which are dashed lines (boson prop- 
agators or interaction lines) such that the diagram separates into two discon- 
nected parts when the dashed line is cut. In the diagram (|4. 1[) the four articu- 
lation lines are highlighted by dots. The diagram (|4.1[) has n a — 4 articulation 
lines. When the articulation lines are cut, the disconnected pieces consist of: 



G 



• cycles, which are closed loops (formed by oriented fermion propagators) 
the vertices of which are connected only to articulation lines. The diagram 
(14. ip has n c = 3 cycles. 

• one (boson) line irreducible parts, which are irreducible in the sense that 
they cannot be separated into two disconnected parts by cutting a dashed 
line. The diagram (|4.1[) has rij = 2 irreducible parts. 

Every diagram which has n a articulation lines, n c cycles and ni irreducible 
parts is such that: 

ni + n c - n a = 1 (4.2) 

This topological property results from the fact that, whenever a boson propa- 
gator is added to a diagram, it is either an articulation line or not. If it is, it 
either adds a cycle or an irreducible part. If it is not, it leaves ni,n c and n a 
unchanged. 

The topological property (14. 2p does not hold for diagrams with open ended 
meson dashed lines, which contain in fact a source factor j at their end point. 
The latter should be counted as an irreducible part and for such diagrams the 
topological relation should read: 

nj + rij + n c — n a = 1 (4-3) 

where nj is the number of open ends (or the number of source points) of the 
diagram. 

In the following we shall refer to one (boson) line irreducible diagrams simply 
as one-line irreducible diagrams. 



5 The particle densities expressed in terms of 
one-line irreducible diagrams. 

In Appendices IA.2I and (|B.2[) the particle densities are expressed in terms of 
connected diagrams T c . In the case of fermions with Coulomb interactions, the 
particle density is given by (IA.22[) : 

d 4 x 1 d i x 2 Y7 — n — r (xi \g\ x) (x \g\ x 2 ) (5.1) 
o {xi \g\x 2 ) 

For nucleons interacting with mesons, the various densities are given by (|B.33[) : 

p a (x) = -tr (x \g\ x)T a + / d 4 x 1 d 4 x 2 -^i — ^ — r ( x i Iffl x ) r « {% \g\ x 2 ) (5.2) 

J b{xi\g\x 2 ) 



In (|5.ip and (J572J) , the connected diagrams are calculated with the unperturbed 
fermion propagator: 

9=7T—T ( 5 - 3 ) 

Ot + ho — (i 
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and r a are the operators denned in (|3.8p and (|3.9p Fl 

From here on we use the formalism applicable to nucleons interacting with 
mesons. Is it however simple to recover expressions applicable to fermions with 
Coulomb interactions. It suffices to set: 

r o = l U a {x) = U(x) p a (x)=n(x) (5.4) 

In the case of nucleons interacting with mesons, the coupling constants are 
included in the r a . As a result, the particle densities p a (x) are multiplied by 
the coupling constants and they are related to the densities n a (x) in the usual 
sense by equations ()B.25|) . In the case of fermions with coulomb interactions, 
the coupling constant e 2 is included in the interaction (|2.5|) . 

We now show that it is also possible to express the particle densities p a (x) 
in terms of one-line irreducible diagrams <E> (U) as follows: 

= -tr (x \G\ x)T a + f d 4 x lC ? 4 x 2 ■ (a* \G\ x) T a (x \G\ x 2 ) (5.5) 

J d(xi\G\x 2 } 



In the expression (|5.5|) , G is the "dressed" fermion propagator: 

o T + ho - [i + u a r a 

and U a (x) is a local potential which satisfies the equation: 

U a (x) = J d A y (x \K ab \ y) {Tr In G" 1 + $ ([/)) (5.7) 

For nucleons interacting with mesons, the boson propagators K a b are defined 
in (|3.12p . For fermions with Coulomb interactions, K a b = K, where K is (mi- 
nus) the Coulomb potential (12. 5|) . The oriented lines of the one-line irreducible 
diagrams $ (U) are the dressed propagators G. When ho is time independent, 
as assumed in this work, the potentials U a (x) and the densities p a (x) are also 
time independent because (x \K a b\ y) depends only on x^ — y^. 

In this formulation, an approximation to the potential U a (x) and the particle 
density p a (x) are obtained by choosing a subset $ (U) of one-line irreducible 
diagrams. Any finite or infinite subset can be chosen. For example, <3? (U) could 
be limited to the following one- line irreducible diagrams: 



*(!/) = 







(5.8) 



J No confusion should arise from the use of V to denote connected diagrams r c and the 
operators 113.81 1. 
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When the set (|5.8[) of one-line irreducible diagrams is chosen, the density (|5.5[) 
is: 



Pa{x) = (x, a) 




+ (x,a) 




(x,a) 




+ (x, a) 



(5.9) 

The representation of the density in terms of diagrams with a slash labeled (x, a) 
is explained in Appendices IA.2I and IB. 21 

The potential (|5 . T[) can be written as follows: 



U a (x) = / d 4 y (x\K ab \y) 



tr (y\G\y)T b - J d 4 Xl d 4 x 2 tr 



5<S>{U) 
S(x 1 \G\x 2 ) 



so that the diagram representation of the potential U a (x) is: 



( Xl \G\y)T b (y\G\x 2 ) 
(5.10) 



U a {x) = (x, a) 



(x, a) — 





+ (x, a) 



+ (x, a) 




+ (x, a) 




I \ + (x,a). 




(5.11) 

In the diagrams (|5.9j) and (|5. 1 1|) the oriented lines are the dressed fermion 
propagators (|5.6| . 

The dressed propagator G can be expressed in terms of the unperturbed 
propagator g by the series: 



G = g- gUTg + gUTgUTg 



(5.12) 



Therefore each oriented fermion propagator G in the diagrams (|5.9p and (|5 . 11 [) 
generates an infinite set of diagrams composed of oriented lines g with any 
number of U insertions. Furthermore, iterations of equation (|5.10[) will generate, 
for each one of these insertions, an infinite set of diagrams connected by an 
articulation line. It is easy to check that the diagrams (|5.1ip generate diagrams 
with a tree structure and that they yield a particle density (|5.2p and a potential 
(|5.7p in which the contribution of each diagram r c is included once and only 
once. 
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In many cases, the ground state densities may vanish because of self-consistent 
symmetries, discussed in Ref.^S , such as spherical symmetry, time-reversal sym- 
metry, isospin or flavor symmetry for example. Self-consistent symmetries de- 
pend on the system which is being calculated, whether, for example, it is in a 
rotating frame, whether it is exposed to an external field, whether the number 
of nucleons is even or odd, whether there is an excess of neutrons or protons, 
and so forth. In this work we do not assume such symmetries. 

Note that the potentials U a (x) and the densities p a (x) are related the simple 
analytic expressions: 

U a ( X ) = -J A ( x \ K ab\ V) Pb (y) Pa (x) = - J d 4 y (x {kJ | y) U b (y) 

(5.13) 

This is in sharp contrast with the potentials which are expressed in terms of 
Kohn-Sham orbits (see Section [T0|) . 

In Section[7]we express the partition function in terms of one-line irreducible 
diagrams and in Section [5] we show that it is a stationary functional of both the 
particle densities p and the potentials U. 



6 Iteration procedure to calculate the particle 
densities and potentials. 

The time-independent single particle hamiltonian hQ+U a T a can be diagonalized: 

(ho + U a T a ) | A) = e\ |A) (6.1) 

The eigenstates |A) are fermion orbits and we distinguish particle orbits \p) 
which have energies e p > \i and hole orbits \h) with energies e/, < p. The 
eigenvalue problem (|6.1[) involves the local and static potentials U a (x) . This 
makes the computation of the fermion orbits |A) considerably simpler and faster 
than, for example, the determination of the poles and residues of the single 
particle Green's function, which involves a non-local and energy dependent mass 
operator. This is the main interest of constructing density functionals. 

In the zero temperature limit, the dressed propagator (15. 6j) can be expressed 
in terms of the particle and hole orbits as in (|C.4|) and (|C.5|) : 

( Xl \G\x 2 ) = (nnlGlrzrz) 
= 9 (n - T 2 )Y,e- {ep -^ (ri - T2) (rib) (P|f 2 )-0(T 2 -rO^VC**-/*)^-™) (n\h) (h\r 2 ) 

p h 

(rrl \G\ rr 2 ) = -J^ (n\h) (h\ f 2 ) (6.2) 

h 

In this representation, the first term of (|5.5|) (but not the second) acquires the 
form of a density of an uncorrelated Slater determinant: 

-tr (x\G\x)T a = J2(h\r)T a (f\h) (6.3) 

h 
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The equations (|5.7[) or (|5. 10[) can be solved by iteration. We begin by choos- 
ing a set $ (U) of one-line irreducible diagrams. Then: 

1. We make an initial guess at the potentials U a (x). Alternatively, we can 
make an initial guess at the particle densities p a (x) and deduce the initial 
potentials from ()5.13[) . 

2. With the potentials U a (x) we calculate the fermion orbits (|6.1[) by diago- 
nalizing ho + U a T a . This yields the dressed fermion propagator G defined 
in fQp . 

3. We express the set $ (U) of one-line irreducible diagrams in terms of the 
fermion propagator G and we recalculate the potentials U a (x) using the 
equations $tT$ or ([5TU|) . 

4. We return to step 2 and we continue the process until convergence is 
achieved. 

The densities p a {x) can then be deduced from (|5.13[) . The equation (|5.2p can 
be used to replace the iteration procedure above by one which yields successive 
approximations to the densities. 

The fact that it is equally easy to calculate either the densities or the po- 
tentials is due to the simple equation (|5.13|) which relates the two. This is why 
we call the eigenstates (|6.1[) ho + U a T a the optimal orbits for constructing sta- 
tionary density functionals. In Section [TU] we show that the use of Kohn-Sham 
orbits is more complicated. 

7 The partition function expressed in terms of 
one-line irreducible diagrams. 

In this section we express the partition function in terms of one-line irreducible 
diagrams. The partition function acquires the form of a stationary functional, 
namely (|7.15p , of either the potentials U a {x) or the particle densities p a (x) . 

In Appendices [Al and [Bl the partition function Z = e w is expressed as a sum 
of connected diagrams T c , calculated with unperturbed fermion propagators g, 
defined in (|5.3|) : 

W = Trlng- 1 +T C (7.4) 

If we wish to express the partition function in terms of one-line irreducible dia- 
grams $ (U), calculated with the dressed propagators G defined in (|5.6p . we need 
to correct for the fact that the one- line irreducible diagrams $ (U) overcount the 
diagrams T c . Indeed, each diagram r c , which can be decomposed into rii (F c ) 
one-line irreducible parts, is included ni (T c ) times in the corresponding one- 
line irreducible diagram $ (U). Consider, for example, the one- line irreducible 
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diagram: 



*(£/) = 




(7.5) 



which is calculated with dressed propagators G. The connected diagram: 





(7.6) 

which is calculated with bare propagators g and which contains nj = 2 one-line 
irreducible parts, is included twice in the diagram $ (U). Indeed either of the 
two irreducible parts of the diagram (|7.6p can be considered as dressing the 
propagator of the other. More generally, a set <!> (U) of one-line irreducible dia- 
grams, calculated dressed propagators G, generates a set of connected diagrams 
r c composed of nj (T c ) irreducible parts and it overcounts them m (T c ) times. 
Its contribution is therefore equal to: 



(7.7) 



There are two other ways to generate connected diagrams. The connected 
diagrams T c generated by the one- line irreducible diagrams $ (U) can also be 
generated by the following set of diagrams, consisting of a fermion closed loop 
interacting once, twice,... with the potential U a T a : 






(7.8) 

In these diagrams the boson lines appended to the closed oriented fermion line 
are interactions with the potential U a T a - The contribution of the diagrams (|T.8[) 
is: 



tr I d A x ( x 



gYaUa ~ ^gr a U a gr b U b + ^gT a U a gT b U b gT c U c - 

= Tr\n(l + gT a U a ) = TrlngQ- 1 (7.9) 

where we used the fact that gG^ 1 = 1 + gT a U a . However, not only does the 
expression (|7.9p include all the diagrams generated by the one-line irreducible 
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diagrams $ (17), but it overcounts them. Indeed, each diagram T c , which can be 
decomposed into n c (T c ) fermion loops, occurs n c (T c ) times in the contribution 
(|7.9[) . which is therefore equal to: 

TrlngG- 1 =J2 n c{T c ) T c (7.10) 

A third way to include the diagrams generated by the set $ (U) is to calculate 
the expression: 





\ J d A x J d 4 y p a (x) (x\K ab \y)p b (y) = i / d A x J d i yU a (x) {x\K ab 1 \y)U b {y) 

(7.11) 

where p a (x) is the density (|5.5p . Not only does the expression (|7.11|) include 
all the diagrams generated by the one- line irreducible diagrams $ (U), but each 
diagram T c which contains n a (T c ) articulation lines, is included n a (T c ) times 
in the contribution (|7.11| . so that: 

\j d A x J d% U a (x) (x \K-i\y) U b (y) = J^n a (T c ) T c (7.12) 

We can combine the results ([777)1 . (fTTTOj) and (f77T2)) to write: 

$ (17) + Tr In gG' 1 -±Jd A xJ d 4 y U a (x) (x \K^\ y) U b (y) 

= K (T c ) + «c (r c ) - n a (r c )] r c = r c (7.13) 

where we used the topological relation (|4.2[) . The expression (I7.13| therefore 
sums the diagrams r c generated by the set $ (U) of one-line irreducible diagrams 
correctly (meaning that each one is counted once). 

The partition function Z = e w can therefore be expressed in terms of one- 
line irreducible diagrams as follows: 

\nZ = W = trlng- 1 + <S> (U) + Tr In gG' 1 - -UK^U (7.14) 

which reduces to: 

W = TrlnG^ 1 + $ (U) - ^UK~ 1 U (7.15) 

In the zero temperature limit, the energy E of the system composed of N 
particles is related to W by the equation: 

W = -j3 (E - fj,N) (7.16) 
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8 Stationary properties of the density functional. 



Consider a variation U a ^ U a + 5U a of the potential U a . The corresponding 
variation of W is: 



SW = 



l 4 x 



6Ua {x) sujx) [ TrlnG_1 + $ (^] - 5Ua 0=) I d 'y ( x \ K ab\v) u b (y) 



(8.17) 

The equation S y W ^ = therefore leads to the equation (|5.7p . In other words, 
the equation (|5.7[) which determines the potential U a (x) is equivalent to the 
equation which states that the functional (|7. 15[) should be stationary with re- 
spect to variations of the potential U a (x) . 

Since U = —Kp, the generating functional (|7.15|) may equally well be re- 
garded a functional of the densities p a - And since ^ = ^7x7 = —K the 
functional (|7.15|) is also a stationary functional of the particle densities p a . 

If all one-line irreducible diagrams $ are neglected, the scheme above reduces 
to the Hartree approximation. If all one-line irreducible diagrams are included 
(a goal never achieved) the scheme yields an exact energy and exact particle 
densities. An approximation is defined in terms of a selected subset of one- 
line irreducible diagrams. In all cases, the fermion orbits are eigenstates of a 
single-particle hamiltonian involving a local potential. 



9 One-line irreducible diagrams identified in the 
path integral. 

The expression (|7.15p has been derived from an analysis of the topology of 
connected diagrams. It can also be derived by an algebraic manipulation of the 
path integral (|3. 15[) of the partition function [3]. Indeed, let us make the following 
change of the integration variable: 

S a (x) -> S' a (x) = S a (x) + U a (x) (9.1) 

where U a (x) is an as yet unspecified local field. The partition function (|3 . 15(1 
becomes, after dropping the primes: 

e W _ J e Trln(d-r + ho-v,+rS+rU+S(®S) + U(<S>U) + U(<S>S)+S(W))-f d l x [^SR- 1 S+SR- 1 U+^U K- X u) 

(9.2) 

Let us define: 

G- 1 =d T + h Q - fi + TU + U($U) (9.3) 

so that: 

Tr In (d T + h - p + TS + TU + S (<I>S) + U ($17) + U (<f>S) + S ($[/)) 

= Tr In G _1 + tr In (1 + GTS + GU ($5) + GS (93)) (9.4) 
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The partition function becomes: 

e W _ e ~Tr\^G-\VK- 1 U J £) ^ g Tr ln(l + GrS+G(7(*S)+GS(*5))-/ d 4 x S+SK' 1 ^ 

(9.5) 

so that: 

W = Tr\nG~ 1 -^UK^U (9.6) 
+ ln f D (S) e^Mi+Grs+GC/^sHGS^S))-/^ (^sk'^+sk" 1 ^ ^ g ^ 



If we compare this to the expression (|7.15[) . we obtain a path integral represen- 
tation of one- line irreducible diagrams: 



$ (U) = In J D (S) 



Tr\n{l+GrS+GU(<i>S) + GS('S>S))-J d l x (^SR- 1 S+SR-^U) 



(9.8) 

To check that the path integral generates only one-line irreducible diagrams, we 
expand Tr In (1 + GTS + GU (<£S) + GS in powers of {TS + U + S ($5)) 

and we derive the corresponding Feynman diagrams as done in Appendices [A] 
and [Bj However, the difference with the diagram expansion of the expression 
(|3.15p is due to the occurrence of the term SK^U in (|9.8p . which generates 
diagrams with open ended interaction lines which are attached to fermion prop- 
agators at one end and to the potential U at the other. These diagrams are 
one-line reducible. The reader can check for himself that, when the poten- 
tial U is chosen to be U = —Kp, these diagrams cancel the one-line reducible 
diagrams (without open ended interaction lines) so that the expression (|9.8p 
generates only one-line irreducible diagrams. 



10 The use of Kohn-Sham orbits. 

We now compare the theory expressed in Sections [5] and [7] to the construction 
of stationary density functionals using Legendre transforms. For simplicity, we 
limit the discussion to systems of fermions with Coulomb interactions, which 
involve only one single-particle density n (x) and for which we can set: 

T a = 1 U a {x) = U (x) Pa (x) = n (x) (10.9) 

Let us add to ho a source term J (x) coupled to the particle density n (x) . A 
stationary density functional is then obtained from the Legendre transform: 

r (n) = W + Jn (10.10) 

The Lagrange multiplier J (x) acts as a local potential. The equation |^ = J 
states that, for a given J, the density is given by n = —^jy and the equation 
4^ = states that J — 0, in which case the density n is the equilibrium or ground 
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state density. The inversion of the equation n = — ^j- consists in finding the 
local potential J (x) which yields a given density n (x) , thereby making T (n) 
a stationary functional of the density. A frequently used way to invert the 
equation n = — ^j- consists in separating the source term J into two parts: 

J=J + Jint (10.11) 

where Jo is chosen such that the exact density n (x) is given by the uncorrelated 
Slater determinant, the orbits of which are eigenstates of the hamiltonian Iiq + Jq 
and which are called Kohn-Sham orbits (see [5],[B], [7] , [8] , [5] , [TU] and references 
therein). The "interacting" part J int is then calculated in terms of the density 
n and the potential J by considering, for example, an expansion in powers of 
the coupling constant, expressed by connected Feynman or Goldstone diagrams, 
calculated with fermion propagators G — g + / 1Q 1 _ M+ j which can be constructed 
from the Kohn-Sham orbits. The interacting part Jj„t can then be expressed in 
terms of one-line irreducible diagrams with, in addition, extra diagrams involv- 
ing for example the inverse density-density correlation function (|C.7j) . Finally, 
when Ji n t is determined, Jo is determined by the equation Jo = — Jint which 
states that, in the ground state, J = 0. 

In section [5] we expressed the density directly in terms of one-line irreducible 
diagrams, which are calculated with a self-consistent potential U (x), determined 
by the equation (|5.7|) . The potential U (x) is then simply and analytically related 
to the density n (x) by the equation (|5.13p . We now show that if we split the 
self-consistent potential U (x) into two terms: 

U{x) = Uo(x) + U int (x) (10.12) 

so that the dressed fermion propagator becomes: 

G _1 = Gq 1 + U int G = G -G U int G (10.13) 

with: 

Go = „ -t 1 ^ TT (10.14) 

Or + h - M + Uq 

and if we choose the potential Uq (x) such that the density n (x) is equal to: 

n(x) = -tr (x\G \x) (10.15) 

then Uq becomes identical to the source term Jo used in the Legendre transform 
method. 

There is no real need to separate the potential into the two terms (|10.12j) . In 
fact we shall see that it makes the theory and the calculations somewhat more 
complicated. We do so however, to compare the present theory to the frequently 
used Kohn-Sham orbits in the Legendre transform method. 

We can express equation (|10. 15|) in terms of the eigenstates of Hq + Uq: 

(ho + Uq) I A) = e\ | A) (10.16) 



16 



and "particle" and "hole" orbits \p) and \h) are the eigenstates |A) belonging to 
eigenvalues respectively e p > fi and e/, < /i. Since Go can be expressed in the 
form (|C.4|) and (IC5|) . the density (|10.15p can be written in the form: 

n(x) = -tr (x\G \x) = ^ (h\ r) {r\ h) 




(10.17) 

The equations (|10.15[) and (|10. 1T|) state the density of the correlated system 
is equal to the density of a Slater determinant composed of the eigenstates of 
+ Uq. That is why the eigenstates (|10.16l) are usually called Kohn-Sham 
orbits. 

10.1 Equations for the potentials U and U int . 

The equation (|5.5[) which expresses the density n (x) in terms of one-line irre- 
ducible diagrams can be written in the form: 

n (x) = - tr (x |G | x) - tr (x \G Q U mt G\ x) - (10.18) 

oil (x) 

where (|10.13p was used. With the choice (|10.15|) of Uq this equation reduces to: 

5<5>(U) 



tr (x\G Q U mt G\x} = 



5U (x) 



(10.19) 

u=u +u tnt 



This is an integral equation for Ui nt which can be solved by iteration, given 
the potential U - Once U int is thus determined, U can be obtained from the 
equation (|5.13[) , which reads: 



U (x) + U int (x) = - J d^x' (x \K\ x 1 ) n (x') (10.20) 
Since K is (minus) the Coulomb potential (|2.5p . we obtain: 

U (r) = / d 3 r' - ,! 2 - s n (f ) - U mt (r) (10.21) 



47r \r — ? 

When U = Uq + Uint the functional 1 5[) becomes: 

W (U) = Tr In G^+Tr In (1 + G U mt )+<5> (U + U int )~ {U + U mt ) K' 1 (U + U mt ) 

(10.22) 

Once U and U int have been determined, we can use (|10.20[) to simplify W (U) 
to: 

W = Tr In Go" 1 - hiKn + Tr In (1 + G U mt ) + $ (U + U mt ) (10.23) 
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10.2 First order approximation. 

Let us first consider the case where the one-line irreducible diagrams $ (U) are 
limited to the first order diagram 



(10.24) 



We consider an expansion in powers of e 2 , considering Go, defined in (|10.14|) . 
to be of zero order. This makes $i proportional at least to e 2 . Let e 2 Ui be the 
first order approximation to Ui n t- 



Ui 7l t (x) = e 2 Ui (x) 
Then, to first order, the equation (|10.19p can be approximated by: 

5<f>i (U) 



tr{x\G e 2 U 1 G \x) = 



SU (x) 



U=U 



We can use (|C.6|) to reduce this equation to: 



so that: 



dV (f| J D|r-)e 2 C/ 1 (r') = §iM 

dUo (x) 



S$i (U ) 



Ui(f) = d 3 / (rlD- 1 ]?) 



(10.25) 



(10.26) 



(10.27) 



(10.28) 



SU (x) 

The equation p0.21|) then yields the potential Uo (r) which generates the Kohn- 
Sham orbits: 




(10.29) 

The diagrams in (|10.29|) are time-independent Goldstone diagrams in which the 
upgoing lines are particle Kohn-Sham orbits and the downgoing line hole orbits. 
The double line is the inverse density-density correlation function (|C.7|) . The 
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first term is the Hartree potential and the second term is the local potential 
which is generated by the Fock term. 

To order e 2 , the functional (|10.23[) reduces to: 

W x (U) = Tr In Gq 1 - ^nKn + TrG e 2 U x + $1 (U ) (10.30) 
We can use (|10.20p to get: 

W l (U ) = TrlnGg 1 - -nKn + TrG (-Kn- U ) +$! (U ) (10.31) 
However: 

- TrG Kn = nKn - TrG U = nU (10.32) 

so that: 

Wi (U ) — TrlnGo 1 + nU + ^nKn + $i (U ) 




(10.33) 



The first order potential (|10.29p and functional (|10.33|) are identical to the 
first order potential Jo (x) and to the functional To + Ti obtained using a Leg- 
endre transform as given, for example, by Valiev and Fernando[8]. 

10.3 Second order approximation. 

Let consider the second order approximation. We limit the one-line irreducible 
diagrams to the first and second order diagrams: 

$1 + $2 = 







(10.34) 

The first order diagram is at least of order e 2 and the second order diagrams 
at least of order e . The second order approximation to Ui„t can be written in 
the form: 

U mt = e 2 C/i + e 4 U 2 (10.35) 

By substituting G = Go — GoUi nt G in the left hand side of equation ()10. 19(1 it 
becomes: 



tr (x \G UintGo\ x) = tr (x \GoU in tGoU m tG\ x) + 



5<P(U) 



6U (x) 



(10.36) 



U=Uo+U int 
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To order e 2 , the equation (|10.36[) reads: 

tr (x\ Goe^Go \x) = S f ) ^ (10.37) 
oU (x) 

which is the result (|10.26p obtained in first order theory, as expected. 

Let us now expand the functional W (U) up to second order. From (|10.23|) . 
and to order e 4 , we obtain: 

W 2 (U) = Tr In G^ 1 - -nKn + TrG {e 2 Ux + e A U 2 ) ~ ^TrG (e 2 ^) G (e 2 f/i) 

+ $1 (U + e 2 f/i) + $ 2 (U ) (10.38) 
For the third term on the right hand side, use (|10.20p to write: 



e 2 J7i + e A U 2 = - J tfV (x \K\ x') n {x') (10.39) 

TrG (e 2 J7i + e 4 U 2 ) = TrG a (-Kn) = nKn (10.40) 

so that: 

W 2 (U) = Tr In G^ 1 + ^nKn - ^TrG e 2 U 1 G e 2 U 1 

+ $i (U Q + e 2 U x ) + $ 2 (U Q ) (10.41) 
We see that W 2 (U) depends only on U\. The third term of (|10.38p is: 

- iTrG e 2 [/iG e 2 C/i = -hr f d 4 x (x \G e 2 U 1 G \ x) e 2 Ui (x) (10.42) 

Z Z J 

We can use f)10.37|) to obtain: 

- -Woe^Goe 2 ^ = -\ [ d 4 x 5 -M^. e 2 U l (x) (10.43) 
2 2 J oUo[x) 

The before last term is, to order e 4 : 

*i (U + e 2 U x ) = $i (U ) + f d 4 x -i^- e 2 C/i (x) (10.44) 

J oUi (x) 

Finally W 2 becomes: 

W 2 = Tr In G l + ^nKn+<S> 1 ((7 )+$ 2 (^o) + ^ / d*x ^J^^ (x) (10.45) 
Since U\ is given by (|10.28p , the last term is: 





(10.46) 



The functional W 2 is exactly the same as the second order functional obtained 
using a Legendre transform as given, for example, by Valiev and Fernando^. 
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11 Comparison to the use of optimal orbits and 
conclusion. 



The potential (|10.29[) and the density functional (|10.45|) are expressed in terms 
of Kohn-Sham orbits (|10.16p which are eigenstates of ho + Uq. Let us compare 
them to the potential and density functional which are obtained with the optimal 
orbits (|6.ip which are eigenstates of Iiq + U. For a given set of one- line irreducible 
diagrams, the potential U (x) is given by (|5.7[) as illustrated in (|5.11|) . When the 
set of irreducible diagrams is limited to the first order diagram (|10.24[) . the 
potential U is given by: 




(11.47) 



and density n is given by the diagrams illustrated in (|5.9[) . namely: 




(11.48) 



The first lines of the diagrams (|1 1 .47(1 are Feynman diagrams whereas the second 
lines are the corresponding Goldstone diagrams. 

The first order potential (|11.47[) . expressed in terms of optimal orbits, is 
simpler than the first order potential (|10.29p expressed in terms of Kohn-Sham 
orbits, because it does not involve the inverse density-density correlation func- 
tion (|C.6j) represented by the double line. But the density (110. 17[) expressed 
in terms of Kohn-Sham orbits is simpler than the density (|11.48p expressed in 
terms of the optimal orbits. However, when the potential is expressed in terms 
of optimal orbits, the density is related to the potential by the analytic expres- 
sion (|5.7|) . which, in the present case of fermions with Coulomb interactions, is 
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simply: 



(11.49) 



Therefore no extra work is required to calculate the density n (f) once the po- 
tential U (r) has been calculated. Furthermore, the second order functional 
(|10.45|) . expressed in terms of Kohn-Sham orbits, has the extra contribution 
(|10.46p which does not occur in the second order functional expressed in terms 
of optimal orbits. When optimal orbits are used, the simple relation (|5.7p be- 
tween the potential U and the density n avoids having to invert the equation 
n(x) = — sj\x) m ^ e Legendre transform method. Furthermore, any (finite 
or infinite) subset of one-line irreducible diagrams can be chosen. The use of 
Kohn-Sham orbits appears to be an unnecessary complication. 

A Diagram rules for fermions with Coulomb in- 
teractions. 

We summarize the Feynman rules in order to fix the diagram notation. A more 
detailed derivation can be found in textbooks [5].|llj. 

A.l Expansion of the partition function in terms of Feyn- 
man diagrams. 

We add a source term j for the boson field luq so that the euclidean action (|2.7|l 
becomes: 



We define an unperturbed partition function in obvious short hand notation: 



I 3 ( w o) 



Tr In (d T + h Q - fi + iu) ) 




and we write the partition function (|2.6[) in the form: 





(A.2) 



so that: 



Wo (j) = Tr In (d T + h - fi) - -jKj = Trlng' 1 - -jKj 



(A.3) 



where g is the unperturbed fermion propagator: 



1 



(A.4) 



9 = 



d T + ho — fi 
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It is the quantity: 



e -W (j) e W(j) = 



-u K 1 w +juo p Tr In(l+giwo) 



/ D((J )ei u ° K ^o+^o 



= /„7Yln(l+giw ) 



(A.5) 



which we expand in terms of Feynman diagrams. In the expression (|A.5[) () 
signifies an integration over the boson field luq according to the expression above. 
The expansion of Tr In (1 + giuJo) reads: 



Tr In (1 + giujo) = Tr giujQ — ~^Tr giu>ogiu>o + 



oo l ,. 

- / d 4 xi...d 4 x n tr (xi \g\ x 2 ) iuj (x 2 ) (x 2 \g\ x 3 ) iuj (x 3 ) ... (x n \g\ x\) iu> 



(A.6) 

The traces Tr and tr are defined in (|2.8[) and (|2.9p . The term of order n gives 
rise to a closed loop formed by n fermion propagators (x \g \ x') represented by 
oriented lines the end points of which are vertices: 




(A.7) 



An oriented line which reaches a point labelled x and which stems from a point 
labelled x' contributes to a diagram the factor: 



(x\g\x>) 



(A, 



When the expansion (|A.6[) is inserted into the expression (|A.5|) and the expo- 
nential e Trlo l 1 + siu<1 ' is expanded in turn, one is left with expectation values of 
products of boson fields: 



(iu (xi) iLu (x 2 ) iuJa (x 3 ) ...) = 



J D (w ) eh UoK Wia ' (iwo (xi) w {x 2 ) ioj (a*) • 
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The expectation value of a single field is: 



(i(j (x)) 



= i 



. SWojj) 
Sj (x) 



i / (fx 1 (x \K\ x 1 ) j (x') 



(A.9) 



The expectation value (contraction) of a product of two fields is: 



(iu) (x) ioj (x')) = (ioj (x)} (iui (x')) + (ioj (x) iuj (%')) 



C 



(A.10) 



and the connected part is the boson propagator: 



(iujQ (x) iLU (x')) 



S 2 W (j) 



c — 



Sj (x) Sj (x') 



(x\K\x>) 



y 



(A.ll) 



Note that the dashed interaction line yields a factor K which is equal to minus 
the instantaneous Coulomb potential, defined in (|2.5p . Expectation values of 
products of three or more fields can be expressed in terms (|A.9|) and (jA.llj) . 
using Wick's theorem. 

Then e~ Wo ^e w ^ is equal to the sum of all distinct unlabeled diagrams 
formed by closed loops of oriented fermion lines and of dashed interaction lines. 
The contribution of a diagram is obtained by the following rules: 

• Assign a label x = (t, r) to each vertex. 

• Each oriented fermion line contributes a fermion propagator (|A.8[) . The 
fermion propagators are written in the order in which they appear as one 
follows the closed loop. For each closed loop we take a trace tr of the 
product of the fermion propagators forming the closed loop. 

• Each dashed interaction line joining two vertices contributes one of the 
factors (|A.9j) or (jA.llj) . depending on whether it joins one or two vertices. 

• A diagram containing n v vertices and ni closed loops is multiplied by the 
factor (Note that in the absence of the source j, the number of 
vertices is always even.) 

• Integrate over the labels x — (r, r) of the vertices and divide by the sym- 
metry factor, which is equal to the number of permutations of the labels 
which lead to an identical unlabeled diagram. 

The contributions of disconnected diagrams factor and their symmetry fac- 
tors are such that the sum of all diagrams is equal to the exponential of the 
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sum r c of connected diagrams. It follows that the partition function can be 
expressed in terms of connected diagrams T c as follows: 

W(j) = W (j) + T c (A.12) 

where Wo (j) is given by (|A.3[) . In the absence of the sources j, the partition 
function is given by: 

In (rre-^-^A = W = Trlng- 1 +T C (A.13) 

where g is the fermion propagator (|A.4|) . An approximation to the partition can 
be defined by a choice of a subset of connected diagrams r c . 

A. 2 Diagram expansion of the particle density. 

Let us add a local source term U (x) to ho. From (|2.3I and (|2.4p we see 
that, upon a variation U (x) — > U (x) + SU(x), we have SlnZ = SW = 
— J d 4 x (rpi (x) SU (x) ip (x)) so that the particle density is given by: 



(A.14) 

(7=0 



When ho is time independent, as assumed in this work, the particle density is 
also time independent: 

(^(x)ip(x)) = (^(r)ip(r)) n{x)=n{r) (A.15) 



However, for any time independent functional F, we have J d 4 x SU (x) jjjjyj = 
P J d 3 r 5U (r) gjjra ■ In order to avoid cumbersome f3 factors, we shall always 
use functional derivatives Trrr-\- 

In the presence of the source, the fermion propagator (|A.4[) becomes: 

g = ar + h-p + u (A - 16) 

From the expression (|A. 13|) of the partition function, we see that the diagram 
expansion of the particle density (|A.14[) is given by: 



n(x) = (x)tp(x)) = - ^ — (Trlng- 1 



(A.17) 



(7=0 



In (|A.16[) . U is the following operator acting in the Hilbert space of a single 
fermion: 

U = [ d 4 x \x)U(x) (x| (A.18) 



It follows that: 

MjJ^j ( x i If I »a) = - (xi \g\ x) (x \g\ x 2 ) 
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SU (x 



■Tr \ng 1 = Tr g \x) (x\ — tr (x \g\ x) 



The first term of (|A.17|) is the unperturbed particle density: 
n (x) = -tr (x \g\x) = ^(ft|r) (r\h) 



(A.19) 




(A.20) 



where we used (|C5|) . For the second term of (|A. 17|) we note that the potential 
U occurs only in the propagators g of the diagrams T c so that: 



8T C 



SU (x) 



= d x\d xi 



$ (xi \g\x 2 ) 



S(x 1 \g\x 2 ) SU (x) 



= — I d x\d xi 



S (xi \g\x 2 

(A.21) 

The second term is the sum of distinct diagrams obtained from the unlabeled 
connected diagrams T c by inserting a slash, labelled x, onto its oriented fermion 
lines. The particle density can therefore be expressed in terms of the connected 
diagrams thus: 



xi \g\x) (a; |5r| ar 2 > 



n (x) = — tr (x\g\x)+ / d 4 xid 4 X2 



5 (xi \g\x 2 ) 

For example, if the only retained connected diagrams are 



x 1 \g\x)(x\g\x 2 ) (A.22) 




r, 



the particle density becomes: 




n(x 





+ 



(A.23) 




(A.24) 

The diagrams (|A.24[) have a slash, labelled x which does not modify the sign of 
the diagram. 
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B Diagram rules for nucleons interacting with 
mesons. 

B.l Expansion of the partition function in terms of Feyn- 
man diagrams. 

For a system of nucleons interacting with several mesons, we add to ho the 
following source terms for the boson fields: 

jS = J (fx [fc (X) „ (x) + JL (X) „. (X) + ,„ (x) flV,.). + i (X) „, (X)] 

(B.l) 

The euclidean action (|3. 14)) becomes: 

/j- (5) = -Tr ln/3 (<9 T + h - y, + (TS) + S (OS)) + ^SK^S - jS (B.2) 
and we write the partition function (|3.15[) in the form: 



,w(j) _ 



= J D(S)e"^ {S) (B.3) 



The unperturbed partition function is obtained by setting to zero the coupling 
constants so that (TS) = S (OS) = 0. In obvious short hand notation: 

e Wo(J) _ e TrU{dr+h -n) J D(S)e-^ SK ' 1S+jS = e Trln(d T +h a -v) + ±jKj 

so that: 

W (j) = Tr In (d T + h — fi) + ^jKj = Tr lng" 1 + ^jKj (B.5) 
where g is the unperturbed fermion propagator: 

o T + h - (i 

In more explicit form, the expression (|B.5|) reads: 

Wo(j) ^Trlng- 1 

1 . 1 . 1 . 1 . 1 . 1 / 1 \ . 

(B.7) 

The traces Tr and tr are defined in equation (|3.6I) . 
It is the quantity: 

e -W ( 3 ) W(j) = jJMP)£ " e 2 = / Trln(l+fl(rS)+ 5 S(*S)) 

/D(5)e-i Sif - ls +^ \ 

(B.8) 
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which we express in terms of Feynman diagrams. In the expression (|B.8[) . () 
signifies an integration over the boson fields S according to the expression above. 

The expansion of Tr In (1 + g (TS) + gS in powers of (TS) + S ($S) 

reads: 

Trln(l + g(TS)+gS($S)) 
= Trg [(TS) + S (<PS)} - \tt 9 [(TS) + S ($5)] g [(TS) + S + ... 

= E - I d i x 1 d i x 2 ...d i x n tr{( Xl \g\x 2 ) [(TS X2 ) + S X2 (*S M )] 

(x 2 \g[ x 3 ) [(TS X3 ) + S X3 ($S X3 )} ... (x n \g\ Xl ) [(TS Xl ) + S Xl ($S Xl )]} (B.9) 

The term of order n is a trace tr of a product of terms the explicit form of which 
is: 

( Xl \g\x 2 )[(TS X2 ) + S X2 ($S X2 )\ 

.r, — Xa) [TV* (X 2 ) + T^ (X 2 ) + r« (8^ a ) X2 + 7r a (X 2 ) (d fl 7T b ) a 

(B.10) 



d T + ho — \i 

The product can be represented by a diagram of the form: 




(B.11) 



28 



in which the oriented (fermion propagator) lines form a closed loop. Each ori- 
ented fermion propagator contributes one of the following factors: 



■r-2 



{xi \g\ x 2 



Xl 



d T + ho — \i 



x 2 ) g<j{3 



•■''1 



(xi \9\X2)TM 



X2H 



Xl 



xi ) fl-w/?7 M 



7T 



•< : i 



= (11 M^rjrWx! 



i 



<9 r + h - fx 
\ 



\ 



/ 



si 



(jGl |5|X 2 ) $a( M fc) = ( Xl 



X 2 U 



Xl 



1 



2 SabcPTcJft 



(B.12) 



<9 T + 7l0 - A* 

Note that the order of the non-commuting operators g and (rS) or S ($5) is 
determined by the direction of the arrow on the fermion propagator. The point 
at which a meson (dashed) line joins a fermion (oriented) line is called a vertex. 
The proper labelling of the vertices is important. Note also that the black dot 
denotes the absence of a derivative of the pion field. 

When the expansion (|B.9[) is inserted into the expression (|B.8[) and the expo- 
nential e Trln ( 1+9 ( rs - l+9S (*' s )) is expanded in turn, we are left with expectation 
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values of products of boson fields which can be evaluated using Wick's theorem. 
The expectation values of single fields are: 



(Si (x)) = / d 4 y (x\K tJ \y) h (y) = (x\K l3 \j 3 ) 



(B.13) 



In more explicit form: 



-d 2 



y^J ia (y) = (x 



8W 

diva 



d 4 y 



-d 2 



■ mi 



yjjwa (y) = (x 



(9 m tt q (x)) = / d A y 



-d 2 



y^j j™ (y) = (x 



Xfl 



dWp 



= d 4 y (x 



(-d 2 + ml) 



' {IV 



m 



yj jw (y) = (3 



Xll 



The expectation values of products of 

(S a (x) S b (y)) = e~ w ^ 



(B.14) 



two fields are: 



s Ja{ x)s Jb{ yf = (s i (x))(s j (y)H(s i (x)s j (y)) c 

(B.15) 



where the connected part is the boson propagator: 

(S a (x) S b (y)) c = , 6 ] W ° ij) ( ; = (x \K ab \ y) (B.16) 
5]a (x) Sj b (y) 

where K ab is the meson propagator (|3.12[) . We can use (|B.7[) to deduce the 
following explicit forms of the boson propagators: 



boson propagators 
(a(x)a(y)) c = (x 



y) = {x 
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a 



(x) u v (y)) c = ( x 



= XjX- 



and: 



(7r Q (x) TT b (y)) c =Sab\X 



1 



-d 2 + ml 



y) =Sab (x 



((9 M 7T Q (x)) 7T b (j/)) = 8ab(x 



y\ = Sab (x 



Xfi 



{■K a (x) (d„7Tb (y))) c = -<5ah ( X 



1 



-i9 2 + 



y) =-8 ab {x\K^d v 



{{dpTta (x)) (dvTTb (y))) c = Sab ( X 



-d 2 



-tK 



y) = -S ab (x 



Xfl 



J, 



(B.17) 



Expectation values of products of three or more fields can be expressed in terms 
(|B.14p and (IB.17|) . using Wick's theorem. 

Then e" w °^'e w ^' is equal to the sum of all distinct unlabeled diagrams 
formed by closed loops of oriented fermion propagators joined to dashed bo- 
son propagator lines at points called vertices. The contribution of a diagram is 
obtained by the following rules: 



Assign a label <j,uj or it to each dashed boson line . Then assign distinct 
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labels to the vertices. An example of a labeled diagram is: 




(B.18) 

Each oriented fermion line, which reaches a point labeled x\ and which 
stems from a point labeled xi , to which a meson dashed line (or two meson 
lines) is attached, contributes one of the factors (|B.12() . 



Each meson dashed line contributes one of the factors (|B.17|) or (|B. 14|) . 
depending on whether it joins one or two vertices. For example, the two 
pion lines in the labelled diagram (IB.18j) . joining the vertices bearing the 
labels X2,X8 and x$, contribute the following factor (remember that the 
black dot indicates the absence of a derivative of the pion field): 



(-)<Ws ( x 2 



1 



-d 2 



2 ^Ms 



X 8 ) (-)K-a 5 { X 2 



<:),< 



1 



-d 2 



-2 ^5 



(B.19) 



The factors are written in the order in which they appear as one follows 
each closed fermion loop and for each loop we take a trace tr of the product 
of fermion propagators forming the loop. 

A diagram containing n v vertices and ni closed loops is multiplied by the 
factor (Note that in the absence of the source terms j a , the 

number of vertices is always even.) 

Integrate over the space-time labels x — (r, f) of the vertices and divide 
by the symmetry factor, which is equal to the number of permutations of 
the labels which lead to an identical unlabeled diagram. 
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The contributions of disconnected diagrams factor and their symmetry fac- 
tors are such that the sum of all diagrams is equal to the exponential of the 
sum F c of connected diagrams. It follows that the partition function can be 
expressed in terms of connected diagrams T c as follows: 

W (j) = Wo (j) + T c (B.20) 

where Wo (j) is given by (|B.7|) . In the absence of the sources j, the partition 
function is given by: 

In (Tre-P^-^A = W = Trlng- 1 +T C (B.21) 

where g is the fermion propagator (|B.6p . An approximation to the partition can 
be defined by a choice of a subset of connected diagrams. 



B.2 Diagram expansion of the particle densities. 

Let us add the following local source terms to ho'- 

r a u a (x) ee rW[/„ (x) + rpu„ (x) + rftu^ (x) + ^ b) u aW (x) (B.22) 

From (|3 . 1 5|) we see that upon a variation U a (x) — > U a {x) + SU a {x), we have: 

SlnZ = 5W = - J d*x (N^ (x)T a SU a {x)N(x)) (B.23) 

so that: 

Pa (x) ee <JVt (x) T a N (x)) = -^^y (B.24) 



The operators r a arc defined in (|3.8p . They include the coupling constants. 
Therefore the particle densities n (x) , in the usual sense of the word, are related 
to the densities p (x) as follows: 

TV (x) EE (N (x) N (x)) = —p a (x) 

n (x) = <7V t (x) N (x)> = — p^o (x) 

n z (x) = (A* (x) a t N (x)> = — p^ (x) 

n^ a (x) = (A* (x) /3r a 757 M iV (x)) = ^-p^a (x) 

9a 

n a (jib) ( x ) = (N (x) EabcTcTTaluN (x)) = 4f%p a {pb) (%) (B.25) 

(In the case of fermions with Coulomb interactions, only one density occurs with 
T a = 1 so that n (x) = p (x).) 
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In the presence of the source term, the fermion propagator (|B.6|) becomes: 

9 = d T + h -p + T a U a (R26) 

From the expression (|B.21[) of the partition function, we see that the particle 
densities (|B.24[) can be expressed in terms of connected diagrams as follows: 



Pa (x) = (Ni {x)T a N(x)) = - —-- (Tr In <T x +r c 



(B.27) 



(7=0 



In (|B.26|) . T a U a is the following operator acting in the Hilbert space of a single 
fermion: 

T a U a = J d 4 x \x)T a U a (x)(x\ (B.28) 

It follows that: 

(xi \g\ x 2 ) = - (xi \g\ x) Y a (x \g\ x 2 ) 



SU a (x) 
S 



Trlng- 1 = Tr g \x) T a (x\ = tr (x\g\x)T a (B.29) 

SUa [X) 

The first term of (|B.27p is the unperturbed particle density: 

p<°> (x) = -tr (x \g\ x) T a (B.30) 

We can use (|C.5|) to get: 

p<°> (x) =J2(h\r)r a {r\ h) = p(°> (r) (B.31) 
h 

For the second term of l|B.27|) we note that the potential U a occurs only in the 
propagators g of the connected diagrams r c so that: 

<^ r c f ,4 .4 ST c s(x 1 \g\x 2 ) f 4 4 sr c 

— — - = d xid x 2 j-, — n — r „ r , r = - / d x x d x 2 — — n — - {x x \g\ x) T a {x \g\ x 2 ) 
oU a (x) J 0{xi\g\x 2 ) 0U a (x) J d{xi\g\x 2 ) 

(B.32) 

The particle density can therefore be expressed in terms of the connected dia- 
grams thus: 

p a (x) = -tr (x \g\ x) T a + / d i x 1 d 4 x 2 — — — - (x x \g\ x) T a (x \g\ x 2 ) 

J b{xi\g\x 2 ) 

(B.33) 

The second term is the sum of distinct diagrams obtained from the unlabeled 
connected diagrams T c by inserting a slash onto its oriented fermion lines, ex- 
actly as in (IA.24j) . except that the slash now bears the label (x,a). An oriented 



34 



fermion line, stemming from a vertex labelled x 2 , reaching a vertex labelled x\, 
and bearing a slash labelled (x, a) contributes a factor equal to: 



= (xi \g\ x) T a (x \g\ x 2 ) 

2 (B.34) 
The slash does not modify the sign of the diagram. 

C Particle and hole orbits and density-density 
correlation functions. 

When ho is time- independent, the unperturbed fermion propagator g^ 1 — (d T + ho — fi) 
is diagonal in the |wA) representation, where the states |A) are eigenstates of ho- 

ho |A> = e A | A) (C.l) 

and where \w) are the (euclidean) time plane wave states: 

(r|u,) = -^e' w d T \oj)=iu>\cj) (C.2) 

so that: 

g^ 1 \u\) — (d T + ho — n) \lo\) = (iuj + e\ — fi) |wA) (C3) 
In the zero temperature limit (3 — > oo, the matrix elements of g are equal to: 

(xi \g\ x 2 ) = (n?i \g\ T 2 f 2 ) = y2 (nri\u\) - — — (wX\ r 2 f 2 ) 

= 9{r 1 -T 2 )Y / e- {ep -^ ){Tl - T2) (n\p) {p\r 2 )~e{T 2 -T 1 )Y j e~ [eh -^ ){Tl - T2) (n\ h) (h\f 2 ) 

p h 

(C.4) 

where the "particle" and "hole" orbits \p) and \h) are the eigenstates |A) be- 
longing to energies respectively e p > /i and < \i. The equal time propagator 
is: 

(rri \g\ rf 2 ) = - £ <fi | h) { h\ ?*) (C.5) 
h 

We shall encounter the so-called density-density correlation function: 

(fi |D| r 2 ) ee -tr / dr (nfi \g\ Tr) (rr\g\ T X r 2 ) = -tr / dr x (t^ \g\ Tr) (rr\g\ T X r 2 ) 
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E — «Pl f 2> (5| h) (h\ fi) (ri|p) + 5) (5b) (p| 5) (?i\h)) 




(C.6) 



The diagram (|G6j) is a time-independent Goldstone diagram in which upgoing 
propagators are particle orbits and downgoing propagators hole orbits. The 
inverse correlation function will be represented in diagrams by a double line: 
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